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1. Introduction 

T-odd effects in hard QCD processes have been attracting our attentions for more than 
30 years, but no experimental verification of the predictions [1-9] has been presented yet. 
T-odd observables change sign under the operation of reversing both the spatial momenta 
and the spins of the all the particles without interchanging initial and final states; see 
Refs. [4, 10] for details.^ In T-invariant theories like perturbative QCD, the T-odd effects 
arise due to the re-scattering phase, or the absorptive part of the amplitudes, which appears 
in the loop level. Such T-odd quantities in hard processes can be predicted in perturbative 
QCD, and should be tested experimentally. 

Since de Riijula et al. proposed to measure T-odd effects as an experimental test of the 
non-abelian nature of QCD in e~^e~ — > T — > ggg with a longitudinally polarized beam [1], 
several theoretical studies have been performed for the quark and gluon processes with an 
electroweak current. They can be classified into three types: 

(i) Three jets in e+e~ annihilation with a longitudinally polarized beam, e^e^ — > qqg [2, 
7,9]. 

^T-odd efTects are sometimes referred to as na'iVe- T-odd [11] or Tjv-odd [9] in order to distinguish them 
from the genuine time-reversal operation T, which exchanges the initial and the final states. 
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(ii) Semi-inclusive deep-inelastic neutrino [3] or longitudinally polarized electron [4, 12] 
scattering, i'hX. 

(iii) Drell-Yan-type process, pp ^* /W / Z + ]et + X . References [5,13,14] considered 
single-spin asymmetry in the Drell-Yan process with longitudinally polarized proton 
beam, while T-odd effects without spin measurement were studied in VF-jet [6, 11] 
and .^-jet [8] events at hadron colliders. 

The absorptive parts of the relevant one-loop amplitudes in these three processes are related 
to each other through crossing [15]. In addition to above three processes, there also exists 
another T-odd observable, the normal polariation in top-quark pair-production at e^e~ 
and hadron colliders [16-21]. 

Observations of T-odd effects in hard processes are a challenging task since they do 
not appear at the tree level. So far, no experimental test has been made for the above 
processes [22-24], even though large non-perturbative T-odd effects have been observed in 
hadron spin physics [23, 25]. We may note that the possibility to observe the perturbative 
T-odd effects in 1^-jet events at the Tevatron run II has recently been pointed out in [11]. 

In this article, we propose a new measurement of the T-odd effects in radiative top- 
quark decays. We study T-odd angular distributions of Pi^-decay leptons in the radiative 
top-quark decay into a bottom quark, a W boson, and a gluon: 

t^h + W+ + g; W+^e+ + i^e. (1.1) 

Due to the large mass, nit = 175 GeV, top-quark decay is not affected by hadronization, 
and hence it can be dictated by perturbative QCD. Even though the correction up to 
0{al) to the total decay width of the top quark is known [26], the correction to the W- 
decay lepton angular distribution in the top quark decay has been calculated only up to 
0{as) [27]. We calculate the absorptive part of the amplitudes for the t — > bWg process in 
the one-loop order 0{a'^), which gives the leading contribution to the T-odd asymmetries. 
The predictions may be tested at future colliders such as the Large Hadron Collider (LHC) 
and the International Linear Collider (ILC). 

The article is organized as follows. In Sec. |2|, we present the lepton decay distribution 
using the density matrices of the t bWg decay and the W ^ decay, and give the gen- 
eral kinematics relevant to our analysis. In Sec. ^, after showing the T-even lepton angular 
distributions, we discuss the T-odd distributions in detail, and study their observability 
at future experiments. In Sec. ^, we consider radiative decays of polarized top-quarks and 
discuss another T-odd observable, the angular correlation between the top-quark spin and 
the decay plane. Section ^ is devoted to a summary. In appendix |A|, we give the absorp- 
tive part of the t bWg decay amplitudes in the one-loop order by using the Feynman 
parameter integral calculation. In appendix we present our results in terms of the loop 
scalar functions. 
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2. t —>■ bW^g decay density matrix 

The decay rate of the process ( |L1| ) can be expressed in terms of the t — > bWg decay and 
the W ^ iv decay density matrices in the narrow width approximation of the W boson, 



(2.1) 



where Tw is the total decay width of W boson, and A, A' = ±, denote the VF-boson 
hehcity. The 3x3 M^-polarization density matrix for the decay reads 



1 dV^ 3 

d cos d d(p Svr 

with the decay branching fraction = B{W — > ii') and 
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Here, the 3x3 matrices are for A, A' = (+, 0, — ), and the polar and azimuthal angles {6, cp) 
of the charged lepton are defined in the rest frame of the W boson, where the z-axis is 
taken along the W momentum direction in the rest frame of the top quark. The x-axis 
{6 = 7r/2, (p = 0) is in the t — > bWg decay plane as explained below. 



Before we show the t bWg density matrix dT^^, , we define the kinematical variables 
for the process 



t{pt,crt) b{pb,ab) + VF+(g, A) +g{pg,(Jg), 
where the four-momenta of each particle are defined in the top rest frame as 



p'^ = (mt, 0,0,0), 



(2.4) 



Pb 



{Eb,Pb sm6,0,pb cos( 

{Eg,Pg,x,0,Pg,z)- 



P'o 



(2.5) 




Helicities of each particle, cr^, cJb, A and a. 



^' Figure 1: Schematic view of the coordinate sys- 
are also defined in the top rest frame. The ^ ^ ^^+^ ^^^^^^ ^^11^^^^ 

z-axis is taken along the W boson momen- ^/-+ _^ decay. 

tum, and y-axis is along q xpb, the normal 

of the decay plane; see Fig. |^. 

We define the dimensionless variables as 



2pfPb 2pfq 2pfPg 



mt 



/ 2Eb 2E\Y 2Eg 



\mt'' mt 



mt 



(2.6) 
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(a) (6) (c) (d) (e) (/) 

Figure 2: Feynman diagrams for the t bWg decay [30]. Tlie top two are tlie tree level diagrams, 
and the bottom six are the one-loop level diagrams contributing to the absorptive part of the 
amplitudes. 



These are the energy fraction of b, W and g, respectively, and satisfy the energy conserva- 
tion condition, zi + Z2 + zs = 2. The kinematically allowed region is given in the (zi, Z2) 
plane by 

2y < 21 < 1 - + y^ 2x < Z2 <l + x^ -y^, 

{z\ - ^y^){zl - 4x2) - [2 + 2x2 + 2y2 _ 2zx - 2z2 + 21^2]^ > 0, (2.7) 

with X = mw /fnt and y = mi,/mt. 

The mass of the 5-quark is kept to be finite (m;, = 4 GeV) for the tree-level calculation. 
However, as we will see later, the effect of the mass is negligible. Thus, for the calculation 
of the T-odd distributions, we take the nib = limit, which simplifies the framework of 
the one-loop calculation. In the case that we ignore the 6-quark mass, there appears a 
kinematical singularity in the Z2 — > 1 + x^ limit, when the 6-quark and gluon momenta 
are coUinear. An infra-red (IR) singularity also exists at Z3 — > 0, where the emitted gluon 
becomes soft. 

Let us now present the density matrix for the t — > bWg decay, dV^^y in Eq. ( |2.1| ). The 
matrix elements of the t — > bWg decay are expressed as 

iMx = ^^t'Vtb uipb, ab) T^" u{pt, at) e* (g. A) 6^ {pg, a^), (2.8) 

where g and gs are the weak and strong coupling constants, t"" is the SU(3) color matrix, 
and Vfb is the Cabibbo-Kobayashi-Maskawa (CKM) matrix element. The tensor T^" is a 
4x4 matrix in the spinor space. The leading contribution to the real part of T'^" comes 
from the tree diagrams [28,29], while the imaginary part appears first in the one-loop 
diagrams. All the tree and the one-loop diagrams needed in our analysis are shown in 
Fig. We give details of our calculation of T^'^ in the appendices. 

Factorizing the color factor and the coupling constants, we define the reduced density 
matrix Hw as 

Y,MxM*y = 4V27rGFasml.\Vtb\^CF • H^y. (2.9) 
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The summation stands for the sum/average of the spins of the particles except W boson 
and the sum/average of colors. In terms of Hw, the density matrix dT^^y is expressed as 



dzidz2 32^2: 



Finally, combining the top- and VF-decay density matrices in Eqs. ( 2.10| ) and (^I^), the 
decay distribution in Eq. (|2.lD is expressed as 



dz\dz2 d COS dd(t) 256 V 27r3 f-; 

A, A 

= K[Fi{l + cos^ 9) + F2(l - Scos^ 9) + F3 sin 29 cos </> + sin^ 9 cos 20 

+ F5 cos 6* + sin 6* cos </> + F7 sin 6* sin + Fs sin 26* sin + F, sin^ 6* sin 2<j) ] , (2.11) 

where K is the factor in front of the summation symbol in the first line. The nine inde- 
pendent functions Fi_9(zi,Z2) are defined in terms of the reduced density matrices H\\i 
as 

Fi = - -1- Hqo + H. 

F2 = o-^OO) 



F 



3 



2 
1 

2^2 





Fq = 






F7 = 


1 

— -f^o+ — + Ho-) , 


H-Q — Ho-) , 


F8 = 


2^ (-H'+o Ho+ + H-o Ho-) 




F9 = 


{H+- - H-+) . 



Fa = \ {H+- + H-+) 
F5 = H++-H^_, (2.12) 

The terms independent of the azimuthal angle, Fi, F2 and F5, are provided from the diag- 
onal terms of the density matrix, while the azimuthal-angle dependent terms are provided 
from the off-diagonal terms, i.e. the interference between the different polarization states 
of the W boson. The terms Fi through Fg are T-even, and the leading contribution comes 
from the tree diagrams. On the other hand, F7 to Fg are T-odd, and they receive the 
leading contribution from the absorptive part of the one-loop amplitudes through the in- 
terference with the tree amplitudes. Parity transformation changes the sign of (j), thus 
Fy^g^g are not only T-odd but also parity-odd (F-odd). Assuming CP invariance, the lep- 
ton angular distribution for the anti-top-quark decay, i — > bW~g; W~ can be 

obtained by changing the sign of Fy^g^g in Eq. (l2ll] ). 



3. Lepton decay distributions 

In this section, we present numerical results for the T-even and T-odd lepton angular 
distributions in radiative top-quark decays. Note that, in our leading-order analysis, the 
T-even distributions Fi_6 are 0{as), while T-odd distributions Fy^g^g are 0{a1). 
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3.1 T-even distributions 



In Fig. ^, we show a contour plot of the function Fi{zi,Z2), which gives the total rate of 
the t — > bWg decay, dT/dzidz2 = ii' (IGvr/S) Fi, after integrating over the lepton decay 
angles. The kinematical boundary given by Eq. ( |2.7] ) for mt = 175 GeV, mw = 80.4 GeV 
and mil = 4 GeV (m;, = 0) is shown by the thick (thin) dotted line. To avoid the IR region 
near Z2 = ^;2max ~ 1-2, we impose the kx cut, 



4 = 2min(p^,p^)(i_cos0bg) > (20GeV)2, 



(3.1) 



where is the angle between the 6-quark 
and gluon momenta in the top rest frame, 
shown by the dashed line. Furthermore, we 
apply the following cut: 



cos 6bn > —0.9, 



(3.2) 



1.2 



1.1 



shown by the dot-dashed line, in order to 
avoid the configuration where the 6-quark 
and gluon jets are anti-collinear. These two 
cuts enable us to define the top decay plane 
spanned by ph and pg, from which the az- 
imuthal angle (p of the decay lepton is mea- 
sured (see Fig. |l|). 

The decay rate is large in the region 
where Z2 is large, because of the collinear 
singularity in the nif, = limit. As the fig- 
ure shows, the effect of the 6-quark mass is 
small for the decay-rate itself, however the 
kinematical boundary as well as the cuts are 
changed slightly by the mass. 

Next, we define the differential asym- 
metries as 




0.9, 



Figure 3: Contour plot of Fi(zi,Z2) in the 
tree level, zi and Z2 are the energy fraction 
of the bottom quark and the W boson, respec- 
tively. The dotted line denotes the kinematical 
boundary; the dashed and dot-dashed lines are 
for the kinematical cuts for kx > 20 GeV and 
cosdiig > —0.9, respectively. The thick contours 
are obtained for mi, — 4 GeV, whereas the thin 
contours are for nib = 0. 



Ai{z2) 



j dziFi{zi,Z2) I j dziFi{zi,Z2) 



(3.3) 



for i = 2 to 9. In Fig. we show the Z2 distributions of the T-even asymmetries A2^— ,6 at 
the tree level for the three zi regions: 2;imm < zi < 0.4, 0.4 < zi < 0.55 and 0.55 < zi < 
-2^1 max) with the same kinematical cuts as in Fig. |3|. The Z2 distributions of Fi for the same 
zi regions are also shown as a reference. In all the figures, predictions for mi, = 4 GeV and 
the massless 6-quark limit are shown by thick and thin lines, respectively. Except for Fi, 
the lines for mi, = 4 GeV and those for = are almost degenerated. Small difference 
in Fi at large Z2 and small zi arises because of the difference in the kinematical boundary, 
as shown in Fig. 0. 
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0.9 1 1.1 1.2 0.9 1 1.1 1.2 0.9 1 1.1 1.2 




0.9 1 1.1 1.2 0.9 1 1.1 1.2 0.9 1 1.1 1.2 



Figure 4: The Z2 distributions of the T-even asymmetries A2 to Aq at the tree level. Three cases 
for the different zi regions with the same kinematical cuts as Fig. || are shown. Thick hnes are for 
rrih — 4 GeV, and thin lines are for m;, = 0. The distributions of Fi integrated for zi are also shown 
as a reference. 

The asymmetries in the polar angular distribution ^2,5 are predicted to be large, more 
than the azimuthal angular asymmetries ^43^4^6- When the VF-boson energy (i.e., Z2) is 
large, the kinematics of the t — > bWg three-body decays becomes close to that of the 
t — > bW two-body decays. Near Z2 = ^;2max; this leads to the well known results: (i) The 
asymmetry A2, which dictates the fraction of the decay to the longitudinally polarized 
W bosons, reaches 0.7. (ii) The fraction to the left-handed W bosons is ~ 0.3, and the 
fraction to the right-handed W bosons is negligible. This corresponds to the asymmetry 

A^ oc — H ~ —H The difference of the factor 2 comes from the normalization in 

Eq. ( |2.12| ). (iii) The ^3,4,6 asymmetries vanish in the large Z2 region, since the interference 
between the different helicity states of the W boson is very small. 

On the other hand, the smaller Z2 becomes, the larger the gluon contribution becomes. 
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Figure 5: The same as Fig. ^, but for tlic T-odd asymmetries ^7,8,9 at tlie one-loop level. 

Due to this gluon contribution, the decay to the right-handed W boson is allowed, even 
in the nib = limit, which causes the deviation from the values in the two-body decay 
process. 

3.2 T-odd distributions 

Let us now turn to the T-odd asymmetries, the main subject of this article. As mentioned 
above, the leading contribution to the T-odd effects in the top-quark decay (0) comes 
from the interference between the tree diagrams and the absorptive part of the six one-loop 
diagrams in Fig. ^. The one-loop amplitudes are calculated in the = limit, however 
the kinematical boundary as well as the cuts are given for nib = 4 GeV. We set the QCD 
coupling constant as Og = as{kTrain = '^0 GeV) = 0.15. The details of our calculation of the 
one-loop amplitudes are given in the appendices. 

Figure |5| shows the asymmetry distributions as Fig. ^ but for ^7.3,9 in Eq. ( |3.3| ). We 
found that the asymmetry A7 is positive at a few percent level, and tends to be larger with 
increasing zi and decreasing Z2- As is also positive but less than 1% in magnitude, and 
is large for the intermediate values of zi and Z2- Ag is the smallest in magnitude and is 
order permill. It takes positive value for large zi and small Z2, but changes the sign by 
decreasing zi and increasing Z2- The dips which appear in the figure are caused by the 
kinematical cuts given in Eqs. ( |3.1| ) and ( |3.2[) . 



In Fig. ^, we show the contribution to the Af asymmetry for 0.55 < zi < zimax from 
the individual one-loop diagrams of Fig. Q in the Feynman gauge. The sum of the diagrams 
(c) and (/), which have the gluon three-point-vertex, gives negative contribution to Aj. 
On the other hand, all the other diagrams give positive contribution to the asymmetry. 
The diagrams (a) and (d) with s-channel 6-quark exchange diagrams give the dominant 
contribution, which make the total asymmetry positive. The diagrams (b) and (e), which 
contain the u-channel 6-quark exchange in the final-state rescattering, are found to give 
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small contribution. On the other hand, the main contribution for comes from the 
diagrams (c) + (/), while for Ag, the contributions from (a) + {d) and (c) + (/) are comparable 
in size. 

3.3 Up-down asymmetry for the LHC experiment 

In order to help finding an evidence of the T-odd asym- 
metries in experiments, we discuss a simple observable for 
the T-odd asymmetry. We define the up-down asymmetry 
A\jD with respect to the top decay plane as 

AuD = [iV(0 <cP<tt)-N{7t <cP< 2n)] /iV,,^. (3.4) 



-0.02 



-0.04 



I I I I I I I I I I I I I I I I _ 
(a)+(d) ^ 

7 




0.9 



1 1.1 1.2 



It is defined as the asymmetry between the number of 
events having the charged lepton momentum with posi- 
tive and negative y component. Aub reflects the property 
of Aj, since sin0sin(/> is positive for < (p < ir while neg- 
ative for TT < (j) <2'K. 

We estimate Aud, and its statistical errors for 820,000 
top-quark signal events which is expected at the LHC one- 
year run with L = 10 fb^^ after the event selection for 
the single lepton plus jets channel — > — > bbWW — > 
bb{iij){jj) [31]. Taking into account the fraction^ of t — > 
bWg events that satisfy the kinematical cuts in Eqs. ( |3.l| ) 
and (|3.2| ), a sample of about 72,000 events for t — > bWg 
followed by W ^ iu would be expected. In Fig. ^ (left), 
we display the distribution of the event sample in the zi-Z2 
plane. In order to see the T-odd asymmetries effectively, we divide the kinematical region 
into eight bins using the jet-energy ordering and the opening angle between the two jets in 
the top rest frame as 



Figure 6: (Left) The contri- 
bution to the A7 asymmetry 
for 0.55 < zi < ziniax from 
the individual one-loop diagrams. 
(a)+(d), (6) + (e) and (c)+(/) 
contributions in Feynman gauge 
are plotted in dashed, dotted and 
dotted-dashed line. Total asym- 
metry is also plotted in solid line, 
as a reference. 



cos 9bg < —0.5, 
-0.5 < cos 9bg < 0, 
< cos 9bg < 0.5, 



(I) zi > Z3 cos Obg < -0.5, (V) zi < Z3 

(II) zi > Z3 -0.5 < cos 9bg < 0, (VI) zi < Z3 

(III) zi > z:i < cos 9bg < 0.5, (VII) zi < z^ 

(IV) zi>Z3 0.5 < cos 9bg, (VIII) zi<Z3 0.5<cos%. (3.5) 

In the figure, the number of events in each bin are given in an unit of thousands. As in 

Fig. ^, a large number of events is expected for the region where both zi and Z2 are large, 

namely (III) and (IV). 

The top and middle plots in Fig. |^ (right) show the up-down asymmetries with expected 

statistical error-bars for each of the eight bins, for the LHC one-year run. The error is 

estimated from 6A = [1 — A'^)/Nsum for each bin. The magnitude of the asymmetry is 

larger for the (I)-(IV) bins than for the (V)-(VIII) bins, and increases with the opening 

^For the total decay width of the top quark, we use the calculation including the 0{as) QCD correc- 
tions [26]. 
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Figure 7: (Left) Estimation of tlie event yields for the LHC one-year run is shown in each bin 



defined in ( |3.5| ). (Right) Up-down asymmetries ^ud defined in Eq. (3.4) for the eight bins (top and 
middle) and Aud for the case without 6-tagging (bottom), cos 6bg is the opening angle between the 
two jets in the top rest frame. Error bars are estimated for the expected event yields shown in the 
left figure. 



angle 6hg, as is expected from the zi and Z2 dependences of in Fig. ^. The asymmetry 
reaches 3% at the bin-(I) where, however, the event yield is not high. 

In the bottom plot in Fig. ^ (right), we also consider the case where the top-pair 
productions are identified without a 6-jet-tagging. In this case, instead of defining y-axis 
by the direction qxpf,,we define the y-axis along q x p^^, where pj-^ is the momentum of 
the jet whose energy is large than the other in the top-quark rest frame. This asymmetry 
corresponds to ^UD for zi > z^ (top) minus ^UD for zi < z^ (middle). Because of the 
cancellation, the magnitude of the asymmetry decreases, but it remains finite even without 
6-jet identifications. 

4. Polarized top-quark decays 

Although we have considered the decay of unpolarized top-quarks so far, the top-quarks 
produced singly by the electroweak interactions at hadron colliders or the top-quark pairs 
produced in e^e^ colliders can be highly polarized. Therefore, it may be useful to analyze 
the polarized top-quark decay. 

In this section, we show that, when a top-quark is polarized, i) there exists another 
type of T-odd observable, the angular correlation between the top-spin direction and the 
top decay plane, and ii) the lepton angular distributions discussed in the previous section 
are modified^. 



First, we discuss another type of T-odd observable in radiative decays of the polarized 
top-quarks, namely, the angular correlation between the top-quark spin and the decay 
plane. 

^We thank the referee of this article for pointing out the existence of another T-odd observable in the 
polarized top-quark decay, and suggesting its relation to the normal polarization in the top-quark pair- 
production. 
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We define the angles between the top- 
quark spin direction and the decay plane 
in the top-quark rest-frame as shown in 
Fig. ^. The z-axis is chosen along the 
T^- momentum direction, and the x-axis is 
chosen along the pb direction in the {pb,Pg) 
plane. The polar and azimuthal angles, Op 
and (pp, respectively, define the direction 
of the top-quark spin st- 

The decay distribution is now charac- 
terized by the two angles as well as zi and 

Z2- 




Figure 8: Schematic view of the coordinate sys- 
tem for the t — > bW~^g decay, where the mo- 
mentum direction in the top-quark rest-frame is 
chosen along the z-axis, with the top-quark's spin 

Sf 



dV 



dzidz2d cos 



GFasm'^x'lVtbl'Cp 
64^/27r3 



Fpi + Fp2 cos 6p + Fp^ sinOp cos (j)p + Fp^ sinOp sin (j)p 



■ (4.1) 



The structure functions i*pi_P4(zi, Z2) are obtained from the t bWg matrix elements 



M-ati which are defined in Eq. (|2.8|), but we now retain the top-quark helicity at instead 
of the W-helicity (A): 



^pi = ^E(l-^+l' + 1-^-1') 



Fp2 = Ie{\M- 



l-M-l^) 



Fp3 = Ie {MIM^ + , 



(4.2) 



The summation stands for the sum of the spins of all the particles but the top-quark, and 
the sum/average of colors. The spin-independent term Fpi is identical to Fi in Eq. ( 2.12| ), 
including the normalization factor. Fpi is T-even and P-even, while Fp2 and Fp-^ are T- 
even and P-odd. Pp4 is T-odd and P-even. The leading-order contribution to the functions 
Fpi to Pp3 comes from the tree-level amplitudes. On the other hand, the leading-order 
contribution to Pp4 comes from the interference between the tree amplitudes and the ab- 
sorptive part of the one-loop amplitudes, just the same as P7,8,9 in Eq. (|1|). Note that 
Pp4 is proportional to the expectation value of the triple product of the three vectors 
{sfijyi Pb), just like P7 is proportional to (sV • q x pb). The corresponding distribution for 
the anti-top-quark decay can be obtained by reversing the sign of Fp2 and Pp3 in Eq. ( [4.1J ) , 
when the CP is a good symmetry. 



We define the ratios of the correlation functions Fi for i 
independent term Fpi as 



A{z2) = j dziFi{zi,Z2) I j dziFpi{zi,Z2) 



P2-P4 to the spin- 



(4.3) 
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Figure 9: The Z2 distributions of the angular correlation functions defined in Eq. (4.3), Ap2.pz at 
the tree level and at the one-loop level, where the three zi regions and the kinematical cuts 
are the same as in Fig. ^. 



Each correlation function corresponds to the expectation value of the component of the 
top-quark spin-vector as 

{st) = ^ (^P3, ^P4, Ap2) . (4.4) 

In Fig. the Z2 distributions of ^P2,P3 at the tree level and Ap^ at the one-loop level 
are shown, where the three zi regions and the kinematical cuts are the same as those in 
Fig. ^. The T-even P-odd asymmetries Ap2 and ylp3 are as large as a few times 10% in 
magnitude, while the T-odd P-even asymmetry Ap^ is less than 1%. This means that the 
top-quark spin lies almost in the decay plane, or, the decay plane tends to contain the top- 
quark spin. The zi dependence of ^P4 is similar to the T-odd lepton angular asymmetry 
Aj in Fig. |. 

Next, we consider the T-odd lepton angular asymmetry Aj again, but in the decay of 
polarized top-quarks. Since the degree of the normal polarization to the decay plane is quite 
small as shown in Fig. |9|, for simplicity, the case that the top-quark spin lies in the decay 
plane, ^p = 0°, is considered. In Fig. |lO|, we show the A-j asymmetry for 0.55 < zi < ^^imax; 
where the spin direction of the top-quark is set at 0p = 0° and 180°. The asymmetry is 
enhanced when Op = 0°, but reduced when 6p = 180°. It follows from the fact that the 
decay amplitude to the right-handed VF-boson is larger for 9p = 0° than for 9p = 180°. 



Finally, we briefly mention T-odd effects induced by the absorptive part of the top-pair 
production amplitudes, which produce the normal polarization with respect to the scatter- 
ing plane. The one-loop calculations have been done for e~^e~ and hadron colliders [16-21], 
however, the degree of polarization is estimated to be quite small. 
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We examine if the up-down asymmetry with respect 
to the decay plane of the top-quarks, studied in this pa- 
per, can contribute to the T-odd asymmetry about the 
scattering plane in the top-pair production process, when 
the production and decay processes are considered in total. 
When the top-quark has normal polarization with respect 
to the scattering plane, because the charged lepton prefers 
to be emitted in the direction of the top-quark spin, the 
expectation value of the inner product of the top-quark 
spin direction and the lepton direction {st ■ pi) is positive. 
On the other hand, considering the T-odd effects in the 
top-decay process, since the asymmetry in Fig. ^ is 
slightly positive, the expectation value of the triple prod- 
uct {st • qx ph) is slightly positive. In addition, since Aj 
in Fig. ^ is positive, the expectation value of {p£ - qx ph) is 
also positive. Therefore, the T-odd effect in the top decay 
process gives positive correction to {st -pi), i.e. additive to 
the original asymmetry due to the T-odd polarization of 
the top-quark normal to the scattering plane. However the 
size should be negligible, because the degree of the normal 
polarization and the T-odd correlation Ap4 are estimated 
to be very small. 

Similarly, we find that the T-odd effect in the top-quark production process provides 
additive but negligible contribution to the T-odd asymmetry in the decay process with 
respect to the top decay plane. 

5. Summary 

In this article, we studied the top quark decay into a bottom quark and a W boson 
accompanied by one gluon emission, and calculated the absorptive part of the t bWg 
decay amplitudes at the one-loop level. We then estimated the leading-order contribution 
to the T-odd asymmetries of the lepton angular distribution in the t — > bWg decay followed 
by leptonic decay of the W boson. 

For completeness, we also discussed the T-even asymmetries at the tree level 0{as), 
and found that the fraction to the right-handed W boson increases in the small VK-boson 
energy region. As for the T-odd asymmetries, the largest asymmetry is predicted for at 
a few percent level, and the other asymmetries {Ag and Ag) are found to be less than 1%. 

We proposed a simple observable ^uDj the up-down asymmetry with respect to the 
top decay plane, which is proportional to Aj. The ^UD asymmetry is predicted to be at a 
few percent level, which may be confirmed at the LHC with 10 fb~^. 

Before closing let us mention that, for the polarized top-quark decays, there exists 
another T-odd observable, the angular correlation between the top-quark spin direction 
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Figure 10: T-odd asymmetry 
Af for 0.55 < zi < zimax in 
the polarized top-quark decays. 
The direction of the polarization 
is parameterized from the W- 
momentum direction. Op — 0° 
(dashed) and 180° (dotted). As 
a reference, the unpolarized case 
is also plotted in solid line. 
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and the top decay plane. However, the size of the T-odd correlation is less than 1%, which 
may be difficult to measure at future colliders. 
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A. t — i> bW~^g decay amplitudes 

In this appendix, we outline our calculation of the amplitude for the t bWg process. 
Note that we present the formalism in the nib = limit, because we performed the one-loop 
calculation only in this limit. The extension to the massive 6-quark case will be given only 
for the tree-level calculation. 

First, we expand the tensor T'^" in Eq. ( |2.8| ) as 

Tf"" = Y,aiTl'° (A.l) 

i 

with the 20 basis tensors; 



rpfJ.a 






= rp^p./mi 


rpfJ.a 

-'ri 


= 9'"'P+/ 


rrit, 


^R6 


= rpUP-^lml 




= j^r4P-/mi 


rpfia 






= -f^rP+/mt, 


^R7 


= p^p?P+/ml 


L3 




rplia 


= PtPbiP-/mt^ 




= ph''4P 




^RS 


= p>tptP+/ml 


rpfia 






= PhPtiP-/T^t-, 




= p'^r4P 




rpfia 


= p^pfP+/ml 


rp^ia 

^L5 




rpHa 


= PbPb4P-/mt, 


rpfJ.Ce 

^R5 


= i^pUp 




rpfia 


= PbPhP+l^t 



(A.2) 



where the chiral-projection operators are P± = 2(1±75). The summation runs for i = {Ll- 
L10,i?l-i?10}. The coefficients Oj are calculated perturbatively, 

cii = bi + iagCi + • • • , (A. 3) 

where bi is the tree-level contribution, and Cj is the one-loop contribution to the absorptive 
part. 
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The 20 coefficients satisfy the following sum rules, because of the gauge invariance of 
QCD {pgT^''' = 0), 

2(1 - y2)aL2 + Z3aL5 + y2aLQ + 2aij2 = 0, 
2aL2 + 2aij2 - z^ciRb - y2am = 0, 
2aLi - 2aL3 + zsan + y2aLS - ^urs = 0, 
2aL3 + 2aiji + 2(1 - 2/2)0^ + z^aRj + y2aR8 = 0, 
2aLi + 4ai2 + 20^4 - zsaig - 2/2aLio + 2aij4 = 0, 
2aL4 - 2aRi - AaR2 + 2(1 - 1/2)0^4 + z^URg + 2/20/110 = 0, 



(A.4) 



where we defined 1/2 = 1 — Z2+x'^ . In appendices A.l, A. 2 and ^ we present the following 



14 coefficients; i = L1-L4:,L6,L8,L10,R1-R4:,R6,R8,R10. The remaining 6 coefficients; 
i = L5,L7,L9,R5,R7,R9 are then obtained from the above identities, Eq. (A. 4)^. 



Counting the number of the physical amplitudes, only the twelve among the 14 coef- 
ficients are independent [4,32]. Using the Dirac matrix identity [32], the terms with Tj^^Q 
and T^"o ^® removed by the following replacements. 



aLi - 


ail 


0-12 - 


aL2 


aL3 - 


aL3 


aiA - 


aLA 


(ILG - 


aie 


ais - 


CLLS 


(^Rl - 


O-Rl 


aR2 - 


aR2 


aR3 - 


aR3 


aRi - 


aRi 


aR6 - 


aR6 


aR8 - 


aR8 



2(2/3 - ziZ2)aLio + -ziaRio, 

1 1 

-Z2y3aLio - ^ysOiJio, 

^(z| - 2x^)aLio + ^Z2aRio, 

^{(zi - ^^2)22 + 2x'^}aLio + ^( 

«L8 + 2:20^10 + a-Rio, 

1 2 1 

-zix CLIO + ^ysflmo, 

1 2 1 

-zix CLIO - -ysa/Jio, 

1 

^y3aLio, 
1 

1/ 

2(2^1 - ^2)aLio + aRio, 

„2 



X OLIO, (A.5) 
with = 1 — Z3 — x'^. 



*To verify our results, we have calculated all the 20 coefficients independently and checked that these 



satisfy the Eqs. (A.4) 
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A.l Tree- level results 

At the tree level, the amplitude has the contributions from two Feynman diagrams (Fig. |2|), 
TCe = 7"^ \ -I'^P- + 7^P--i -7". (A.6) 



The decomposition in terms of T-^^ in ( |A.2| ) gives 

bhi = bis = -bm = 2xb, bLi = -bie = 2xt, -bL2 = &ii2 = xt + Xb, (A.7) 
where xt = m1/{—2pt • Pg) and xi, = m^/2pfe • pg. 

For the massive 6-quark case, two more components, 

TZ = g^^"P^/mt, Tj:^, = rrP-/m, (A.8) 

with the coefficients 6mi = 2yxh, 6m2 = —y{xt + x^) and y = rrib/mt, must be added to 
Eq. 



A. 2 One-loop results 

At the one-loop level, the absorptive part emerges from the six diagrams for the t — > bWg 
decay, shown in Fig. ^ We write the one-loop coefficients in Eq. (A. 3) as the sum of these 
diagrammatic contributions, 

c^ = ci^+^+^ + ^+c^+^^■ (A.9) 

The analytic expressions of the coefficients are obtained for each diagram by performing 
the standard Feynman integrals. Our expression contains functions with a one-parameter 
integral, which can easily be evaluated. 

In the next appendix, we also show the results of Cj in the loop scalar function method 
as an alternative expression. We checked that the numerical results of the two calculations 
agree completely. 

With the color factor Cp = 4/3, Ca = 3 and Ci = Cp — Ca/2 = —1/6, the one-loop 
coefficients for each diagram are found as below; 



diagram- (a) 



diagram- (6) 



diagram- (c) 



.(a) _ oJa) _ _Ja) _ Ja) _ _r,Ja) _ Cf_ 

2 



2c^2 — ~C_j^3 ~ '^Rl ~ '^^R2 — ~^b- (A. 10) 



-2cS = 4cS = -2cS = cS = 2cS = -4c2) = C,Xb. (A.ll) 



41 = 2cg = -cS = 2cS = = -241 = ^Xb (In + In x,) , (A.12) 



-16- 



where e = mg/mt- The gluon mass nig is introduced to regulate the IR singularity. We 
keep e only in the singular parts and take the e — > limit elsewhere. 
• diagram- (d) 



-LI 



L2 
Cf 



-Cf 



2- V) 

Xb{2 - Z2)ho - (3 - Z2)l21 + ^7^/22 



Xb{l- X^)/io 



(1 - X')l21 



y2 



h2 - x^y2 ihs - -^34) 



-2c' 



id) 
R2 



Cf 



Xbil - x^)ho - (2 - x'^)l2i + 



y2 



I22 



= Cf [hi - ^22 - 2/2 (^32 - ^33)] ■ 



(A.13) 



Here, Imn is defined by the integral 



i [1- 



Z2t + xH'^Y 



(A.14) 



• diagram- (e) 



xi,\n{zfxb) - (1 + zi)JuQ + {z2 - x'^)Jiu + y^2i3 - Z1L2 



» __1 (e) , 



■'L2 



Cl3 = Cl [Xb\n{zfxb) - ho + y2l21 - (1 + 2:1)^110 + {Z2 - X^)Jiii + ?/|Jl21 



2''^i+ 2 



-2/2 



^l(«^121 — ■^122) + "^211 — 2J212 H T^J21^ 



+2/2^211 7, <^212 - Z1L2 



■'L4 



Cl 



JlW - {Z2 - X^) Jill - 2/2(1 - Z2) J12I 



2/2(2/2 + 2a;2 



-Cl 



2x6 + Jiio - (1 + a;^) Jill - ^2(2 + zi- ^2)^121 + 



- J122 



2/2(2 - 2/2) 



J122 



^ 2/2(2/2 + 2x^) ^^^^ _^ £li£l±Mi,2 



^1 -2/2 



2^f2/2 

2;i - y2 



i^3 



-Cl [4Jiii - 22/2 ( J121 + J122 + 2J212 + J213) + 2/i ( J223 + 2 J314)] , 



Clio = <^i2/2 [2J122 - 2/2 (2J133 + J224)] 



-iii 



.(e) 

-Ll 



+ Ci 



'10 



2/2J 



110 



^ (J212 — J213) 



-R2 



I + Z2 — 

-2Jiii + 2/2 J121 + (2 - 22) J211 J212 



-Ciy2 [J211 — J212] , 



-714 



-Cl [Jiio — Jill + 2/2 ( J121 — J122)] , 



^R6 = ^1 [-^llO - 2Jlll + 2/2 (J12I + J212)] , 

cSg = Cl [2Jiio - 2y2 (J121 + 2J211) + 2/i (J222 + 2J313)] , 
cmo = -Ciy2 [4Ji2i - 2J122 - 2/2 (2Ji32 + J223)] , 



(A.15) 
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where Jmni -^n are defined as 



t'^dt 



Lfi 



dt 



[y2t + Zi{l- t)]^ "A 1 - ^2t + XH^ ) ■ 



In 



(A.16) 
(A.17) 



diagram- (/) 



„(/) 



.(/) _ 
'L2 ~ 



.(/) 



■'L4 



'L6 



'LIO 



.(/) 
„(/) 



J/) 



J/) 

-Re 



Ca 
2 



4 



2 



Xb In - X6 - (1 + zsj^iio + o -^iii + n "^m 



^^2^2 



J1' 



122 



, ^2(1 + 2:3(^2 - SZ3) , 
+2/2.^212 



2 -^213 2(^2 - .3) 
-,2, 



2 - Z2 + 2a;^ 



111 



a;t(lne +\nxb) - Xb\n{z^Xb) + Xt + Xb + {I + Z3)J[iq - ^ 

y2{z2 + Z3)j, .27/ 2/2 w .X^y2 J, , 3. . ^r/ w 
^ -^121 + X ^2^122 - y -^211 + ^--^213 + l^\y2 + ^3)^2 " 2/2^3-^^3 

' 1 / 2 N or/ ^^2- 22:3-23;^ , , ^2 (32:2 + 2:3 -4X^) w 

a^fcln (233;&) -^6 -22;3-^iio 2 2/2^121 H 2 1 

n ^ 7' _L ^2(3 - 2^3 - x^) , 

-2/2(1 - 2:3)^211 + 7, -^212 



122 



£3(y2__3£3)^/ 
2(2/2 - ^3) ' 



2/2^3 



2/2 - Z3 



■L' 



Ca 
' 2 



xt (Ine^ + Inxb) + + J[-^q - {z2 - x^) J^i - ^^(^ + ^2) ^/^^ 



y2(l + 2z2--x2) , 2/2(32/2 - ^3) r' 

' 7t "^122 ' Tu \~^2 



Xb 

2 



2(2/2 - 2:3) 
) (in + In Xb) +xt + Xb + J'no — 



■L' 



yhz 
y2 - zz 



1 + Z2 + X' 



•^111 



^2(2^2 + 2:3) 



121 



^2(2^2 -^3) J, 2/2(2-22) J, , 2/2(1 + X^) , , 37/2+223 „ , r/' 

1 ^ •^122 7, -'212 -I 7, -'213 "I 7^ ^2 - 2/223-t^3 



=L8 = ^ [4^ni - 2y2 (3J{22 + 2J212 + ^13) + yl (2^3 + J'22^ + J224 + 2J314)] 



.(/) 



— 2/2 [2J(22 - 2/2 (2J(33 + 424)] , 

,(/) , '^^ 



+ 



,(/) + ^ 



'L2 



'10 



y2 
2 ' 

2/2 + 423 , 



- y i ^3(^^(21 - -^122) + J211 - (2 + 2/2)^12 + J: 



213} 



110 



(2/2 - 22:3) J(n - f J{2i 



2 ^110 



2/2(1 - 22) , x^y2 ( , , x' 

7 '^211 7~ {'^'^212 ~ 'J213) 



Jill ~ 22/2 («^121 ~ "^122 + "^211 ~ "^212)] ! 



.(/) _ '^A 



~2' ["^110 "^ni ^^2 (■^121 <^122)] 5 



Ca 



[j[io - 2-^ni + 2/2 (-^122 + <^2i2)] 
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— —Ca [J'uo — 2^2 (-^121 + ^211) + 2/2 {^133 + ^223 + ^Sis)] 



C 



(/) _ '^A 



RIO 2 



-iy2 [4J(2i - 2J[,, - 2/2 (2J{33 + ^3)] > (A.18) 



where J^„^ and are given by replacing zi to Z3 in Jmni and L„, in Eq. ( A.16|) and 
Eq. ( A.17| ), respectively. 



We note that the sum of the IR singular terms from the diagrams (c) and (/) is 
exactly proportional to the tree-level amplitude, therefore they do not contribute to the 
T-odd distribution. 



B. Loop scalar functions 

As a check of our calculation, we calculate the one-loop coefficients in terms of the loop 
scalar functions, the Passarino and Veltman's B, C, D functions [33]. 

For each diagram, we assign the masses and the momenta of the scalar function, fol- 
lowing the FF notation [34], and take only the imaginary part of the functions. In this 
assignment we explicitly present the 6-quark and gluon mass, mi, g, for clarity, even though 
we take the massless limit in our analysis. 

• diagram- (a) 

Defining Bi = Im Bi{m'^,ml; p^^) for i=0,l with p^^ = {ph + Pg)'^, the coefficients are 
expressed as 

-cS = 2cS = -eg = 41 = -2^1 = Cfx, [Bo + B,] /vr. (B.l) 

• diagram- (6) 

Defining cf* = lmCi{mg,ml,ml; plg,p'^,pl) for i=0, 11, 12, 21-24, the coefficients are 
expressed as 

_Sb) _ 2^b) _ (b) _ (b) _ (b) 

= Ci [ - Co - 2Cn + C12 - C21 + C23 - 2C24/plg] ^t/^, 
cf^ = Ci[-Co- 2Cii + C12 - C21 + C23] ml/TT. (B.2) 

• diagram- (c) 

Defining cj-*^^ = Im Cj (m^ , , ; pig, pi, p^) for z=0, 11, 12, 21-24, the coefficients are 
expressed as 

^Ll ~ ^'^L2 ~ ^L3 ~ ^Rl ~ ^^i?.2 

= Ca[Co - Cn - 2C21 + 2C23 - 12C24M'J ml/ATT, 
c'll = Ca [2Co + 4Cii - 3Ci2 + 2C21 - 2C23] ml /An. (B.3) 

• diagram- (d) 
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Defining C^^'^^ = ImCi{rng,m'f,ml-, Pt,q^,plg) for i=0, 11, 12,21-24, the coefficients are 

expressed as 

eg = -2cg = -CpXbi- (2 - Z2)Co - (3 - Z2)Cu + {Z2 - x')Ci2 - C21 - x2C22 + -22C23 

- 2C24M'] my IT, 
eg = -Cfx^I- (1 - x^){Co + Cii) - x2(C22 - C23) - 2C24/mf] m?/7r, 
eg = -2eg = -CfX6[(1 - x2)Co + (2 - x2)Cii - {z2 - x^)Ci2 + C21 + x2C22 - Z2C23 

+ 2C24/m?] mf /tt. 



c)?3 = -<^^^fe [-C11 + C12 - C21 + C23] m|/7r. (B.4) 



• diagram- (e) 

Defining d\^^ = lmDi{ml,mlml,ml;pf,q'^,pl,pl,plg,{q + pgf) for i=0,ll-13,21- 
27,31-313, the coefficients are expressed as 

eg =Ci [ziDo + (1 + zi)Dn - Du + 2(D27 + ^312 - ^sis)/^?] mf/Tr, 

eg =Cl [ - ZlDo - (1 + Zi)Dn + i?12 + I?21 + (^2 - x2)I)22 - (1 - Z2 + x2)Z)23 " 2L>24 
+ ^1^25 - {Z2 -Z3- 2x^)D2Q - 2D27/ml - X^D^2 - D34 + + Z2D3e 

- zsDsr - (1 - zi - 2x^)D3s + {1 - zi - x^)D^q - {z2 - Z3)D^w 

- 6(1)312 - D3i3)/mJ\ mt/2Tr, 

eg =Ci [ziDo + (1 + 2zi)Du + (1 - - 2^2 + 2x^)0^2 - (1 - Z2 + x^)i2Di3 - D23) 
+ (2 + Zi)D2l + x'^D22 + {l-zi- 3Z2 + 2x'^)D2i - (4 - 2^2 + X^)D2^ 
+ {2z2 - 3a;2)D26 + -D31 - ^2-034 - (1 + ^3)^35 + x^D^q - D^s) + 23^37 

- (1 - zi - x2)D39 + {l-zi + Z2- x2)L>3io + 4(I?27 + -D311 - D3i3)/ml] mf/n, 
('lI =Ci [ - ^11 + {Z2 - x'^)Du + {l-Z2 + a;2)L>i3 - D21 - x^D22 + Z2D24 + £'25 

- {Z2 - x2)L>26 - 2(£>27 + £'3i3)/m?] m^/TT, 

eg =Ci [Dn + (1 - 2z2 + x^)Di2 + {z2 -Z3- x'^)Dis + 2^21 + x^D22 

- (1 - Z2 + X^)D23 - 2Z2D24 + {Zi - Z3)D25 + (1 - Zi)D26 + D35 - Z^D^-j 
+ X^Dss + {l-Zi- x'^)D3g - Z2D3W + 2(^27 - ^312 + 3r>3i3)/mf ] m^/TT, 

eg =Ci [2(L>i2 - Dis + D24) + £'22 + i^23 - £'25 - 3£>26 + £'36 - D38 + £'39 - i^sio] 
X 2mj/TT, 

ego =Ci [-D23 + D26 + D3S - Dsg] 2m^/7r, 

eg =Ci [ - ziDo - {2zi + Z2- x^)Dn - {1 - zi - 2z2 + 2x^)Di2 + (1 - ^2 + x^)Di3 

- D21 - X^D22 + Z2D24 + ZSD25 - (l- Zi- X^)D26 

- 2(21)27 + ^311 - ^313)/"^?] mj/7T, 



='^1 [-21^0 + {2ZI +Z2- X^)Dii + (1 - Zi - 2Z2 + 2x^)Di2 + (2 - Z2)D; 



21 



(1 - Z2 + x2)(Z)i3 - D23) - {Z2 - 2X^){D2A - D26) - (3 - 2^2 + X^)D25 + £'31 
Z2L'34 - (1 + -23)£'35 + x'^iD^e - D^s) + Z3D37 - {1 - Zi - X^)D39 
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+ {l-zi + Z2- x^)D3io + 2(2^27 + - 3L»3i3)/m?] mt/2Tr, 

('m =Ci [-D21 + D24 + D25 - D26] mj/n, 
41 =Ci [^11 - -^12 - D25 + D2e] mt/TT, 
('m =Ci [-Du + 2Di2 - Di3 + D24 - D26] mt/Tv, 

=Ci [-Du + - 2D21 - D23 + SD25 - Du - + 2L>3io] 2m^/7r, 
41o =Ci [D23 - 2D25 + D26 + Dsg - Dsio] 2mt/7r. (B.5) 
• diagram- (/) 

Defining D^^^ = ImDi{ml,mlml,ml; p'f,q'^,pl,pl,plg,{q + Pbf) for i=0,ll-13,21- 
27,31-313, the coefficients are expressed as 

=-Ca[- 2Z3D0 - 2(1 + 23)^11 - (2 - 3z2 + 2x^)Di2 + (2 + zi - 2z2 + 2x^)Di3 

- 3D21 + {Z2 - 3x^)D22 - (2 - 3Z2)D24 + 3ziD25 + (3 - 2zi - 3Z2 + 3x^)D26 

- X^D32 - Du + Z2D36 -{I-Z3- X^)Ds8 + ziDsio - 2(2^27 + ^3i2)/m2] 

Cl2 = - Ca[2z3Do + 2(1 + Z3)Dn + (3 - 4z2 + 3x^)Di2 - {2 + zi - 2z2 + 2x^)Dis 

+ 5D21 + 3X^D22 - 4Z2D24 - 5ziD25 + 4(1 -Z3- X^)D2Q + QD27/mf\ mi/SiT, 

ClI =-Ca[- 2z^{Do + 2Dii) + (Z3 - ^1)^12 + ziDii - (1 + 2^3)1^21 - 3a;2D22 

+ 2(1 -zi + x^)D2A - 2(1 - 2zi -Z2 + x^)D25 -{I-Z3- x^){3D2e + L'310) 

- A31 + ^2^34 + Z1D35 - x^DsG - 2(7L»27 + 5Dsn)/mj] mj/An, 

ClI =-Ca[- 2Dn + 2(z2 - x^)Di2 + (1 - ^2 + x^){Du + D23) - 2D21 - 2x2l>22 

+ 2z2i^24 + (2 - 3Z2)D25 + (^2 + 2x2)^26 + ^35 " ^1^37 + X^ D38 
+ {1-Z3- x2)L>39 - Z2D310 - 2(2L>27 - 5L>3i3)/m2] m^/47r, 

Cl6 =-Ca [2Dn + (3 - 5z2 + 3x'^)Di2 - 2{zi -Z2 + x'^)D^3 + AD21 - {z2 - 4x'^)D22 

+ 2;i-D23 + (2 - 5^2)i^24 - 2(1 + 2Zi - Z2)D25 - (3 - 2zi - AZ2 + 5x'^)D26 

+ x'^D32 + D34 - D35 - Z2D3Q + Z1D37 + (1 - Z3 - 2x'^)D38 -{I-Z3- x2)£>39 

- {zi - Z2)D3w + 2(3£>27 + £'312 - D3iz)/mf\ mj/AiT, 

Cl8 =-Ca [2(£'l2 - £»13 + -D24 - £'25) + i^22 - £'26 + ^36 - i^31o] r4/lT, 

Clio = -Ca [£'23 - -D26 - £'38 + -D39] m|/7r, 

c^ll =-Ca [2z3Dq + 2(1 + 2z3)L>ii - (2 - 2zi + Z2- 2x'^)Di2 - (2 + zi - 2^2 + 2x'^)Diz 

+ 5L>21 + 3x2£>22 - 4(2;2i524 + ZxD2^) + (1 - 2:3 - x2)(3L>26 + -D310) + £>31 

- 2;2-D34 - ziL»35 + x'^D^Q + 2(5£)27 + D3u)/mf\ rr4/AiT, 

CrI =-Ca[- 2Z3D0 - (11 - 4zi - 5z2 + x^)Du + (2 - 2zi + Z2 - 2x^)Di2 
+ {2 + zi- 2z2 + 2x^)Dri - (5 + ^2)^21 - 5x^1)22 + (5^2 + 2x^)D24 

- (1 - 6zi -Z2 + x^)D25 - 5(1 - zs - x^)D26 - l2D27/mJ\ m^/STT, 

c^g =-CA[-Dn+ D12 - 2{D2i - D24)] m^/27r. 
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Ca[Du - Du + 2(L>25 - ^26)] mt/27r, 

Ca [ - -Dii + 2Di2 - Dis + D2i - D2^] mt/2Tr, 

Ca [ - -Dii + i^ia - 2(L'2i - ^25) - £'34 + £'35] m^/TT, 

Ca [-D2Z + 2D25 - ^26 - £'37 + ^310] mi/'K. 



(B.6) 
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